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We investigate the zero temperature chiral phase transition in an SU{N) 
^ • gauge theory as the number of fermions Nf is varied. We argue that there 

exists a critical number of fermions Nj, above which there is no chiral symme- 
try breaking or confinement, and below which both chiral symmetry breaking 
pH I and confinement set in. We estimate Nj and discuss the nature of the phase 

transition. 



X 

^ ' An SU{N) gauge theory, even at zero temperature, can exist in different 



phases depending on the number of massless fermions Nf in the theory. The phases 
are defined by whether or not chiral symmetry breaking takes place. For QCD 
with two or three light quarks, chiral symmetry breaking and confinement occur 
at roughly the same scale. By contrast, in any SU{N) gauge theory, asymptotic 
freedom (and hence chiral symmetry breaking and confinement) is lost if the number 
of fermions is larger than a certain value (= 11 A/2 for fermions in the fundamental 
representation) . 
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If the number of fermions Nj is reduced to just below llN/2, an infrared 
fixed point will appear, determined by the first two terms in the beta function. By 
taking the large N limit or by continuing to non-integer values of Nf |^], the value 
of the coupling at the fixed point can be made arbitrarily small, making a pertur- 
bative analysis reliable. Such a theory with a perturbative fixed point is a massless 
conformal theory. There is no chiral symmetry breaking and no confinement. 

As Nf is reduced further, chiral symmetry breaking and confinement will set 
in. There have been lattice Monte Carlo studies of the Nf dependence of chiral 
symmetry breaking 1^]. For example, Kogut and Sinclair found that for = 3 
and Nf = 12 there is no chiral symmetry breaking, while Brown et. al. Q have 
found chiral symmetry breaking for = 3 and Nf = 8. In this paper we will 
estimate the critical value Nf at which this transition occurs. We then investigate 
the properties of the phase transition for Nf ~ A^^;. 

Our discussion will parallel an analysis of the chiral phase transition in QED3 
and QCD3 1^, Q. In a large Nf expansion it was found that an appropriate effective 
coupling has an infrared fixed point with strength proportional to 1/A^y, and that 
as Nf is lowered, the value of the fixed point exceeds the critical value necessary 
to produce spontaneous chiral symmetry breaking. It was argued that this critical 
value is large enough to make the 1/A^^ expansion reliable. 

An A^^ dependence similar to the one we describe here has been found in 
A^ = 1 supersymmetric QCD [^. This theory is not asymptotically free for large 
enough A^^-, and has an infrared, conformal fixed point for a range of Nf below a 
certain value. 
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The Lagrangian of an SU{N) gauge theory is: 



where ip is a set of Nf 4-component spinors, the T"" are the generators of SU{N), 
and ^(/i) is the gauge couphng renormahzed at some scale fi. The renormahzation 
group (RG) equation for the running coupling is: 

d 

IJL—ai^) = (3(a) = —ba^ifi) — ca^ifi) — da^{^) — ... , (2) 

Ojl 

where a{^) = g'^{fi)/4Tr. With the Nf fermions in the fundamental representation, 
the first two coefficients are given by 

b = (lliV - 2Nf) (3) 

c = ^ - lONNj - 3^^iV/j • (4) 

The theory is asymptotically free if 6 > {Nf < ^-N). At two loops, the theory 
has an infrared stable, non-trivial fixed point if 6 > and c < 0. In this case the 
fixed point is at 

a* = - - . (5) 
c 

Recall that the coefficients b and c are scheme- independent while the 
higher-order coefficients are scheme-dependent. In fact one can always choose a 
renormalization scheme such that all the higher order coefficients are zero, i.e. they 
can be removed by a redefinition of the coupling (change of renormalization scheme) 
g' = g + Gig'^ + G2g^ + ... . Thus if a zero, a^, of the f3 function exists at two loops, 
it exists to any order in perturbation theory . Of course if the value of is large 
enough, there could be important higher order corrections to the Green's functions 



of physical interest. Indeed, their perturbation expansion might not converge at 
ah. In addition, non-perturbative effects, such as spontaneous chiral symmetry 
breaking, could eliminate even the existence of the fixed point. If the quarks develop 
a dynamical mass, for example, then below this scale only gluons will contribute to 
the /? function, and the perturbative fixed point turns out to be only an approximate 
description, relevant above the chiral symmetry breaking scale. 

For Nf sufficiently close to UN/ 2, the value of the coupling at the infrared 
fixed point can be made arbitrarily small. The RG equation for the running coupling 
can be written as 

Mog f 1 - J- - J- log f " - , (6) 



a q:()u) a* \a{fi) {a — a 
where a = a{q). For a, a{ii) < a* we can introduce a scale defined by 



A = fx exp 

so that 



-1 , /a* — a(fi)\ 1 
log ' 



6 a* V / ba{iJ,) 



i=Mog(f) + i-log(^). (8) 

a V A/ a* Vq;* — a J 

Then for q' S> A the running coupling displays the usual perturbative behavior: 

while for g <C A it approaches the fixed point a*: 

a « "* ■ . (10) 

As Nf is decreased, the infrared fixed point a* increases. We will suggest 
here that the breakdown of perturbation theory, described above, first happens 



4 



due to the spontaneous breaking of chiral symmetry, and that the phase transition 
can be described by an RG improved ladder approximation of the CJT ||^] effective 
potential. It is well known |^] that in vector-like gauge theories the two-loop effective 
potential expressed as a functional of the quark self-energy becomes unstable to 
chiral symmetry breaking when the gauge coupling exceeds a critical valueQ : 

_ _ 2tiN 

where C2{R) is the quadratic Casimir of the the representation R. Thus we would 
expect that when Nf is decreased below the value A'^^ at which a* = Oc, the theory 
undergoes a transition to a phase where chiral symmetry is spontaneously broken. 
The critical value A'^^ is given by 



Nf = N 



\ 25iV2 _ 15 I V ; 



For large N, Nj approaches 4A^, while for = 3, A'^^ is just below 12. Note that 
this is consistent with lattice QCD results |^], which suggest that 8 < Nj < 12. 

Is this simple analysis reliable? After all, it could be that when a* is as large 
as Oc the perturbative expansion for the CJT potential has broken down. To address 
this question we provide a crude estimate of the higher order corrections to the CJT 
potential. An explicit computation of the next-to-leading term (or equivalently the 
next-to-leading term in the gap equation) |jl^ for a* ~ Oc, produces an additional 



factor of approximately e = This is the factor remaining after the appropriate 
renormalizations are absorbed into the definition of the coupling constant. From 



more general definition of tlie critical coupling is that the anomalous dimension of ipip 
becomes 1 [HI. 
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equation (|r^) we see that 

6(1-1^) 

For QCD, e ~ 0.19. If higher orders in the computation produce approximately this 
factor, the perturbative expansion of the C JT potential may be reliable^. The same 
may be true of the various Green's functions encountered in the skeleton expansion 
of the CJT potential. 

We next explore the nature of the chiral phase transition at Nf = N'j and 
its relation to confinement. It is useful to consider first the behavior in the broken 
phase Nf < N'j ( > ac)- Here each quark develops a dynamical mass S(p). For 
Nf — > N'j from below {a^ ac from above), S(p) can be determined by solving 
a linearized Schwinger-Dyson gap equation in ladder approximation. For momenta 
small compared to A, the effective coupling strength is a*, while for momenta above 
A it falls according to equation (^). The resulting solution for S(0) is |15] 



E(0)p.Aexp|^=^l . (14) 

etc , 

The behavior of S(p) as a function of -p will be discussed shortly. 

Once the dynamical mass S(j)) is formed, the fermions decouple below this 
scale, leaving the pure gauge theory behind. One might worry that this would 
invalidate the above gap equation analysis since it relies on the fixed point which 
only exists when the fermions contribute to the /3 function. This is not a problem, 

however, since it can be shown that when S(0) ^ A the dominant momentum 

^It is worth noting that in condensed matter physics one can often (though not always) obtain 
useful information from the Wilson-Fisher expansion in a parameter that is set to one at the end 
of the calculation. 
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ange in the gap equation, leading to the exponential behavior of equation (14), is 



S(0) < p < A. In this range, the fermions are effectively massless and the coupling 
does appear to be approaching an infrared fixed point. Note that the condition 
11(0) ^ A is indeed satisfied for Nf sufficiently close to A^^^;. 

Below the scale S(0) the quarks can be integrated out; thus the effective 
/3 function has no fixed point and the gluons are confined. The confinement scale 
can be estimated by noting that at the quark decoupling scale S(0), the effective 
coupling constant is of order Oc- A simple estimate using equations (|2|)-(^) then 
reveals that the confinement scale is roughly the same order of magnitude as the 
chiral symmetry breaking scale. When Nf is reduced sufficiently below A^^^; so that 
a=K is not close to Oc, both S(0) and the confinement scale become of order A. 
The linear approximation to the gap equation will then no longer be valid, and it 
will probably no longer be the case that higher order contributions to the effective 
potential can be argued to be small. 

It is interesting to compare the behavior of the broken phase for Nf near 
Nf to the walking technicolor gauge theories discussed recently in the literature 
[13|. We have argued here that for Nf just below Nf, the dynamical breaking is 



governed by a linearized ladder gap equation with a coupling just above Oc- As 
the momentum p increases, a{p) stays near a^, (it "walks") until p becomes of order 
A, and only falls above this scale. It can then be seen that the dynamical mass 
S(p) falls like 1/p (i.e. the anomalous dimension of ipil^ is ~ 1) for S(0) < p < A and 
only begins to fall more rapidly ( like ) at larger momenta. This is precisely 
the walking behavior employed in technicolor theories and referred to there as high 
momentum enhancement. In that case, however, there was no IR fixed point to 



keep the (3 function near zero and slow the running of the couphng. It was noted 
instead that the same effect would emerge if the (3 function was small at each order 
by virtue of partial cancelations between fermions and bosons. 

From the smooth behavior of the order parameter 5](0) (equation ([T^)), it 
would naively appear that the chiral phase transition at Nf = Nj (a* = Oc) is 
second order. In this paper we will use the phrase "second order" to refer exclusively 
to a phase transition where the correlation length diverges as the critical point is 
approached from either side. In other words, there is a light excitation coupling 
to the order parameter that becomes massless at the critical point. In the broken 
phase, this mode would be present along with the massless Goldstone modes. In 
the symmetric phase, all these modes would form a light, degenerate multiplet, 



becoming massless at the critical point [11|. 



We examine the correlation length by working in the symmetric phase and 
searching for poles in the (flavor and color-singlet) quark-antiquark scattering am- 
plitude, computed in the same (RG improved, ladder) approximation leading to 



equation (14). The analysis is similar to that carried out for QED3 If the 
transition is second-order, then at least one pole should move to zero momentum as 
we approach the critical point (i.e. the correlation length should diverge). We take 
the incoming (Euclidean) momentum of the initial quark and antiquark to be q/2, 
but keep a non-zero momentum transfer by assigning outgoing momenta q/2 zizp 
for the final quark and antiquark. Any light scalar resonances should make their 
presence known by producing a pole in the scattering amplitude (when continued 
to Minkowski q'^). 

If the Dirac indices of the initial quark and antiquark are A and p, and those 



of the final state quark and antiquark are a and r, then the scattering amplitude 
can be written (for small q) as TxpariPiQ) = ^Xp^arTip, q)/p'^ + where the ... 
indicates pseudoscalar, vector, axial-vector, and tensor components, and we have 
factored out l/p"^ to make T(p,q) dimensionless. We contract Dirac indices so 
that we obtain the Bethe-Salpeter equation for the the scalar s-channel scattering 
amplitude T{p,q), containing only t-channel gluon exchanges. If ^ q^, then q^ 
will simply act as an infrared cutoff in the loop integrations. The Bethe-Salpeter 
equation in the scalar channel for p ^ A is: 

where A is the scale introduced in equation (^). (Note that contributions from the 
integration region k'^ > A^ are suppressed by a factor p^ / h? , and a falling a{k).) 
The first term in equation (jl^) is simply one gluon exchange. We have used Landau 
gauge (^ = 1) where the quark wavefunction renormalization vanishes to lowest 
order. Because of the existence of the fixed point, it is a good approximation to 
have replaced a{p) and a{p — k) by a^, at momentum scales below A. 



For momenta p > q , equation (15) can be converted to a differential equa- 



tion with appropriate boundary conditions. The solutions have the form 

T{p, q) = A{q) y + B{q) V , (16) 



where rj = ^Jl — a^^/ac- The coefficients A and B can be determined by substituting 
this solution back into equation (^). This gives: 

-2.^1 -r,f (^) 



(1+??) 1 _ /.2\r; 



(17) 



and 

^ = {2 ■ (18) 

Note that there is an infrared divergence in the Umit — > in both equations ( [T7| ) 
and (|l|). That this is an infrared divergence rather than a pole corresponding to 
a bound state can be seen from the fact that the divergence exists for arbitrarily 
weak coupling (a* — > 0). In fact, it can already be seen at order a1 in the one-loop 
(two gluon exchange) diagram. As required by the KLN theorem [0], this infrared 
divergence will be cancelled in a physical scattering process by the emission of soft 
quanta. 

If we denote the location of the poles of the functions A and B in the complex 
plane by ^q, we have 

We see that there is no pole that approaches the origin = as a,, — > Oc- Thus 
the correlation length does not diverge, and the transition is not second order (It 
is not conventionally first order either since the order parameter vanishes continu- 
ously at the critical point.). Note that the behavior of the zero temperature chiral 
phase transition is different from the finite temperature case due to the presence 
of long-range gauge forces. At finite temperatures, gluons are screened, and thus 
there are only short-range forces present and only conventional first or second order 
transitions are possible. 

To conclude, we have argued that as the number of quark flavors, Nj, is 
reduced, QCD-like theories in four dimensions undergo a chiral phase transition at a 
critical value Nj (equation (|l2|)). For Nf < Np chiral symmetries are spontaneously 
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broken, while they are unbroken for Nf > Nj^. We have explored the nature of 
the chiral phase transition, arguing that it can be described using the QCD gap 
equation in ladder approximation (equivalently the two-loop approximation to the 
CJT potential). We have also argued that even though the order parameter vanishes 
at the critical point, the correlation length does not diverge (i.e. the phase transition 
is not second order). The critical behavior described here is similar to that found 
in QED3 and QCD3 j^, ^. We have, of course, not proven that higher order 
corrections to our computation are small. Further study of this question as well as 
lattice Monte Carlo studies of the zero temperature phase transition would help to 
confirm or disprove our conclusions. 
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